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Growth of interfaces during vapor deposition are analyzed on a discrete lattice. For a rough surface, 
relation between the roughness exponent a, and corresponding step-step (slope-slope) couplings is 
obtained in (1-1-1) and (2-1-1) dimensions. From the discrete form and the symmetries of the growth 
problem, the step-step couplings can be determined. Thus a can be obtained. The method is applied 
to the linear EW as well as to the non linear growth equation, the Kardar-Parisi-Zhang equation in 
all dimensions. Further, exact exponents for the fourth order linear and nonlinear terms relevant to 
the growth with non conserved noise are obtained in (1+1) and (2+1) dimensions. 



It has been shown recently [1] that the kinetic consid- 
erations on a growing surface can lead to various growth 
terms in a growth equation representing a conserved 
growth from vapor. A growth equation helps in under- 
standing the roughening of a growing interface. It rep- 
resents dynamics of height fluctuations due to the depo- 
sition noise. The non-linear behavior in stochastic equa- 
tions is of general interest. The mode coupling due to 
the non-linearity can lead to the nontrivial physical ef- 
fects that are crucial in understanding many stochastic 
processes [2,3]. However, in the presence of non-linearity, 
absence of analytical solutions in many cases does not 
allow a conclusive understanding of the roughening phe- 
nomenon. In the present work, slope components are 
discrctized (i.e. steps). The discrete form depends upon 
the height difference between the nearest neighbor sites. 
It represents step at a given site and is proportional to 
the local slope. In the presence of non conserved noise the 
height difference fluctuates within the correlation length. 
We refer to these fluctuations as the step fluctuations. 
Coupling conditions between the step fluctuations can be 
determined from the discretization of the relevant terms 
in the growth equation and the symmetry associated with 
the growth problem. These conditions determine the ex- 
act asymptotic behavior of the growth. In particular, in 
a growth, roughness along the interface is measured in 
terms of the exponent a and the time evolution of the 
height-height (h-h) correlation in terms of the exponent 
z. Knowledge of a leads to the determination of z, in the 
case of conserved growth from the scaling behavior of 
noise and in the case of KPZ equation from the Galilean 
invariance [3]. The method is further applied to fourth 
order linear and non linear terms considered to be impor- 
tant in understanding the growth phenomenon. Since we 
obtain exact results the renormalization behaviour of non 
linear terms can be obtained without ambiguity. 

In the following we develop the method for obtaining 
the critical exponents for the second order growth equa- 
tion, both conserved and non conserved. Same method 
is later extended for fourth order terms. 



A linear equation representing interface motion normal 
to the surface can be obtained in the frame of reference 
moving with the interface velocity by considering inter- 
planer hopping of adatoms on the interface with a bias 
for downward or in-plane hopping toward step edge [1]. 
It has the form 

— i^oV^/i + ry (1) 

where, i^q explicitly depends upon F, and 77 is the 
noise due to the randomness in the deposition flux. 
It has the correlation given by < ri{x,t)ri{x' ,t') >— 
2DS{x~ x')(5(i — t'). The angular brackets denote the 
average of the contents. Eq. (1) is known as Edward- 
Wilkinson (EW) equation [4]. The lowest ordered non- 
linear correction to EW equation was introduced by Kar- 
dar , Parisi , and Zhang [5]. The resulting equation, 

dh 

— ^y^-^^h + \{Vhf+'n (2) 

is known as KPZ equation. This is a non-conservative 
equation. 

The KPZ equation has large number of applications in 
growth [2,3]. Its variants are also useful in analyzing var- 
ious physical situations [6]. The directed polymer repre- 
sentation of the KPZ equation [7] is one of the most stud- 
ied equation for understanding the non-linear dynamics 
of the interfaces. Being non-linear an exact solution of 
the KPZ equation is not possible analytically. However, 
many physically relevant quantities are obtainable from 
the analysis of the equation that can be measured ex- 
perimentally. In particular the scaling exponents that 
characterize the growth can be measured experimentally. 
By applying scaling transformations, x — *■ bx, t — *■ b^t, 
and ft, —> 6" ft, to the KPZ equation and demanding that 
the equation remains invariant under these transforma- 
tions one obtains, z = 2, z + a = 2, and z = 2a + d 
for the equation in d dimensions. Since KPZ equation is 
obtainable from the noisy Burger's equation [8,2,3] using 
the transformation v — > Vft, where v is the velocity field. 
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Galilean invariance is implied in the KPZ equation. The 
coefRcient A explicitly appears in the Galilean transfor- 
mation rendering it constancy under renormalization [3]. 
As a result , 



z-\- a. ' 



(3) 



is the relation valid in all dimensions. One can measure 
OL from the height-height (h-h) correlations, 



(4) 



where, correlation length ^(t) ~ t^l^. In the limit 
a; ^ 0, / ^ 1. Time exponent /3, where z = a/ (3 can 
be obtained by measuring the width over a substrate of 



(4)' It 



length L as, W2 = ^ T,xiKx,t) - hf = L'^°'g 

can be shown that [2] for small times 11)2 ^ t^^ ■ 

In 1+1- dimensions a has been determined exactly us- 
ing perturbation expansion [5] and also using Fokker- 
Planck equation with fluctuation - dissipation theorem 
[2]. In 2-1-1- dimensions perturbation theory fails since 
renormalization process indicates existence of only strong 
coupling regime. The dimension dc = 2 is a critical di- 
mension for KPZ equation. For d > 2, renormalization 
reduces A to zero , signaling EW behavior in the weak 
coupling regime. The strong coupling regime has a dif- 
ferent behavior than EW equation. Thus, for d > 2, 
a phase transition is observed for KPZ equation [2,3]. 
To obtain various exponents in dimensions d > 2, non- 
perturbativc methods are employed. There are computer 
simulations [9-11] and theoretical methods [12-14]. Most 
of the theoretical approaches suggest that there is an up- 
per critical dimension beyond which strong coupling van- 
ishes indicating weak coupling EW like behavior. Sim- 
ulations [9] suggest that there is no upper limit on the 
dimension. Using effective large distance field theory sub- 
ject to a few phenomenological constraints, it is shown 
that a—2/h and 2/7 in 2+1 and 3+1 dimensions respec- 
tively. This claim is not supported by large scale com- 
puter studies [10]. It is shown [10] that the values of a are 
not rational numbers as predicted in reference [12]. Us- 
ing a pseudo spectral method [11], numerical solution of 
the KPZ equation in 2+1 dimensions is obtained. Based 
on this method, a is obtained from the saturated width 
(0.37 + 0.02), height-height correlation (0.38 + 0.02) and 
structure factor (0.40 + 0.02). In the present work we ob- 
tain exact values of the exponents of the KPZ equation. 
This is done by discretization of the KPZ equation and 
from the symmetry requirements, evaluating the step- 
step couplings. A step fluctuation at site i is deflned to 
be {hi — hi+i). We will refer to it simply as 'step'. It is 
proportional to the local discrete slope at site i. 



Consider a one dimensional substrate with a lattice 
constant a such that a > Id, where Id is the diffusion 
length. The physical lattice constant will be < a. We 
consider growth on this substrate where atoms are de- 
positing on the physical lattice and follow the relaxation 
rules as depicted by a conservative growth equation. In 
1+1- dimensions such an equation can be written as, 



dh d , 



(5) 



where, J is the local particle current at x. Eq.(5) is 
written in the frame of reference moving with the average 
velocity of the interface defined by the incident flux F. 

In 1+1- dimensions, the (h-h) correlations in the dis- 
crete form are 



G{n) = {{Jh - h,+^f) 



(6) 



We assume that the correlation length ^ is very large com- 
pared to the length na. We define step at i as 

Sxi = hi — hi+i (7) 

The local slope is then —Sxi/a. In terms of steps we have 

G(2) =< Sx^ + dx^^-^ + 2SxiSx,+i > . (8) 

Let < Sx^ >=< 6xfj^-^ >= (5^ and , < 6xiSxi+i >= sS"^ 
where s is the coupling between the steps around i. The 
distribution for Sxi is always symmetric around zero and 
time independent [15] for an ensemble average. 

In the limit ^ — > 00 Eq. (4) reduces to G{x) = ca;^" 
where constant c = G(l). Hence Eq. (8) can be written 
as 



-|2a 



2 + 2s 



(9) 



where G(l) = in the discrctized case. Coupling s 
uniquely determines a. Thus, for s=-l/2,0,and 1, a is 0, 
0.5, and 1 respectively. 

We analyze the growth equations in 1+1- dimensions in 
the discrete form to obtain a for various terms. Consider 
EW equation, 



dh ,^ ^ , 

— = vo{dXi_i - 5xi) + r)i 



(10) 



The first term on the right side represents second deriva- 
tive in terms of local steps. In order that the steady 
growth continues , at any instant the surface must be 
characterized by the fluctuations in the second term 
i.e. various values of second derivative. These must 
be consistent with the requirement that < Sx^ >= 5'^ 
and that corresponding distribution of must be scal- 
able. From the analysis of the linear second order growth 
equation [2,3] a = 0.5. Hence s must be zero for the 
EW equation, s = correspond to Gaussian distribu- 
tion of slopes with average fluctuation of 5"^. The av- 
erage squared difference between consecutive slopes is 
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{5xi — = 2(1 — s)S'^ and the average squared 

sum is ^(^££i±|£i^^'^ = 1(1 + s)d^. For given d^, the 

volume in the parameter space formed with these two 
parameters, squared difference and sum of consecutive 
slopes, decides the number of configurations available 
for the second derivative in Eq. (10). This volume is 
maximum for s = 0. This result can be generalized to 
any higher order term by noting that the steps in the 
present case are proportional to the particle current cor- 
responding to the EW term. Thus by following above 
argument for any other conservative term in 1+1- di- 
mensions , the condition for maximum number of con- 
figurations is same as < JiJi+i >= 0. Thus, when all 
the values —1/2 < s < 1 are scalable in 1-1-1- dimen- 
sions, in the presence of non conserved white noise the 
EW term drives the surface morphology to maximize the 
configurations for the linear second order term. 

For KPZ equation, the scaling transformations lead to 
a = 1/3, i.e. s — —0.206299. Renormalization is ex- 
pected as a result of the non linearity. The KPZ expo- 
nents in 1-1-1- dimensions are obtained without resorting 
to perturbation approach [2]. Here we obtain by refer- 
ring to the step-step coupling. For a scalable surface 
(5^ is constant. For a given 6"^, each scalable value of s 
will correspond to certain distributions of slopes on the 
surface. Any value of s lying between -1/2 and 1 is pos- 
sible under renormalization since Eq. (2) does not indi- 
cate any condition that will prefer any particular value 
of s. Under these conditions, the surface is expected to 
exhibit the maximum range of fluctuations in the local 
slope Sxi consistent with the average of 5"^. Fig. 1 dis- 
plays the distributions of slopes for s = and s = 1 with 
= 1. Distribution for s = 1 is sharp and symmet- 
ric on both sides and for s = it is Gaussian. For non 
zero value of s, the average squared difference and sum of 
consecutive slopes constrain the fluctuations at a given 
site. e.g. for s = 1 the average slope difference is zero 
while sum is 6"^ giving the sharp distribution as in Fig. 
1. For any scalable s 7^ 0, a double peaked distribution 
is obtained. However, the Gaussian distribution provides 
the maximum range of slopes. Hence, for KPZ equation 
the renormalization drives the system toward s = with 
a = 0.5. 

In 2-1-1- dimensions, we need to consider the steps in 
two directions. We define, 5xi_j — h{i, j) — h(i + 
^UiJ = Khf) - Hhj + 1), Sx^^ij = h{i - - h{i,j), 
dyij-i = h{i,j — 1) — ft,(z, j),....over a square lattice. 
We consider isotropic situation so that, < dx^j >=< 
dyf j >= 6^ and the linear coupling constant in the two 
directions is same, s. In 2-|-l- dimensions, apart from 
the linear couplings in the two directions, cross-coupling 
< SxijSuij > also contributes, a can be expressed in 
terms of the cross coupling through the following rela- 
tion. 



< {h{i,j - 1) - h{i + >=< ih{i,j - 1) - h{i,j) f 

+ {hii,j)-hiz + l,j)^ 
+2{h{i,j - 1) - h{i,j)){h{i,j) - h{i -h 1, j)) 

This can be written as, 

G{V2) =< 5xIj_, > + < 5yfj > 
+2 < Sxij-iSyij > 

This leads to 

2" = 2-|-2g (11) 

where, < dxij-i6yij >= qS"^. For —1/2 < q < 0, a 
varies between to 1. Prom Eqs. (9) and (11) we obtain 

S = l + 4q + 2q^ (12) 

Consider a set of N number of pairs of slopes such 
that the average cross coupling over this set is qS'^. One 
can characterize such a set by a distribution of slopes, 
along with a distribution of difference and sum of the 
slopes. All these distributions are needed simultaneously 
to obtain the same coupling q by averaging. Same is 
true in the case of linear coupling. In the analysis to 
be followed we will consider set of pairs of steps while 
averaging over cross or linear couplings. 

To analyze various terms in 2-|-l- dimensions, we con- 
sider a two dimensional substrate with N number of sites 
on a space rcnormalizable lattice. TV is large enough 
for averaging purpose. The growth is assumed to have 
reached steady state with correlations developed over 
large distances. Various configurations are developed on 
the surface in statistically significant number so that the 
spatial averaging can be considered to be equivalent to 
the ensemble averaging . In 2-1-1 dimensions given a ref- 
erence height, morphology can be defined by specifying 
the N pairs of slopes, one at every lattice point. In order 
to specify it in terms of linear couplings, a pair in x- di- 
rection and a pair in y- direction are necessary . Thus a 
quadruplet of slopes at a given site is needed. The mor- 
phology defined by these quadruplets for a scalable sur- 
face is such that the average over linear step coupling is s 
and that over the cross coupling is q consistent with Eq. 
(12). However, the relevant growth term in the growth 
equation demands certain relationship between s and q, 
which helps in obtaining the growth exponents. We ana- 
lyze the EW equation in 2-1-1- dimensions. The EW term 
in discrete form is {dxi-ij —Sxij +dyij-i —6yij). Steps 
in this term form a quadruplet of steps at a given site i. 
Thus there are N quadruplets defining the frozen EW 
surface. This surface represents solution of the Eq. (1). 
This set of steps can be considered to form 2A'' cross cou- 
plings (SxijSyij) and {Sxi^i_jdyi_j^i). The morphology 
of the surface is defined through the specification of set 
of quadruplet defined at every site, i.e. by specifying the 
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value of V^/i at every point. Let the cross coupling con- 
stant be ■ This surface can also be described by linear 
couplings defined over the same set of quadruplets over 
N sites. Let Sew be the linear coupling characterizing this 
set of 2A^ linear couplings {6xij6xi-ij) and (dyijdyij-i) 
over A'' number of sites. Let 2sew/'2'(lew = 1 + P ; where 
— 1 < p so that < 5xi^i_j5xi_j >pf + < Syi^ijSyij >n= 
-(1 +p){< dxijdyij >N + < Sxi-ij6yij-i >n)- The 
subscript N denotes the averaging is over N sites and the 
negative sign ensures the directions in which the slopes 
are measured are consistent. This implies that, the lin- 
ear coupling Sew over the given surface representing so- 
lution of EW equation is same as cross coupling on 
another surface where slopes from the quadruplets de- 
fined over the original surface are scaled by \/l+p while 
keeping the 5^ constant. These two conditions cannot be 
simultaneously satisfied over the same original surface 
for any non zero p. Hence p must be zero. This implies 
Sew = Qew From Eqs. (11), (12) wc get a = for the 
EW term in 2-|-l- dimensions. Thus in 2-|-l- dimensional 
case, the surface morphology is characterized by scalable 
V^/i term for the condition s = q. 

Consider KPZ term in 2-|-l- dimensions. The discrete 
form is {6'^Xij) + (S'^ytj). This form suggests that on the 
N sites there are N pairs of slopes {Sxij , Syi_j). These 
pairs will be characterized by a cross coupling qkpz- This 
set of pairs will have the appropriate rotational symme- 
try. The surface morphology is defined over A'^ sites with 
this set of N pairs. Thus in 2-f 1- dimensions, over N 
sites it is necessary and sufficient to have only N pairs 
of slopes to define the morphology of the KPZ surface 
as compared to 2N in the case of EW surface. Corre- 
sponding to this surface there will be 2N linear couplings 
(Sxi^ijSxij) and (Syij^iSytj) each being defined over 
N sites and characterized by Skpz ■ We will consider set of 
N quadruplets {Sxi-i,j,Sxij,Syij-i,dyij) defined over 
the N sites (same as in the case of EW term) over which 
these linear couplings are defined. The set of cross cou- 
plings is also defined over the same set of quadruplets 
since the symmetry of the problem allows any of the cross 
coupled pairs to be chosen from a given quadruplet. Let 
2skpz/qkpz = 1 + p. As is discussed in the case of EW 
equation, we identify linear couplings 2skpz same as cross 
coupling q on a surface where steps are scaled by a factor 
of y^T^^ and 5^ remains unchanged. This is possible 
only for p — 0. Hence, Skpz = qkpz/2. Using this relation, 
Eqs. (11) and (12) we obtain 

a = 0.35702 (13) 

for this surface, accurate up to 5 digits. One can compute 
it with arbitrary accuracy. Note that for KPZ equation, 
the surface can also be defined using 2A'' pairs of slopes as 
in the case of EW surface. However, for KPZ surface this 
2N pairs of slopes are sufficient to define a morphology 
but not necessary. This will lead to the well known result 



of EW surface for KPZ case in the weak coupling (in this 
case small A in Eq. (2)) limit. 

In 3-1-1- dimensions, we consider a cubic lattice. For 
the additional degree of freedom we define Szi jj^ = 
h{i,j, k) — h{i,j, 1) as the step in the z direction. The 
discrete form of the KPZ term is (S^Xij^k) + {^'^Vi,3,k) + 
(J^Zijjt). There is a triplet of slopes at every site on 
the surface. It is required that average cross coupling 
over N sites for any two directions must result in to 
same q value, qkpz- We can choose a set of A'' pairs of 
cross coupling between Sxij^k and Syij^k characterized 
by qkpz- The slope Szij^k must be chosen so that its 
cross coupling with Sxij^k or 5yij,k is qupzS"^ - To sat- 
isfy this condition we consider at every site the average 
value \{5xij^k -|- Syij^k) as one of the slopes in the cross 
coupling. The other slope is Szij_k chosen in such a way 
that corresponding coupling will result in to a set of cross 
couplings that averages over A'' sites to qkpzS^- This en- 
sures the requirement that cross coupling between slopes 
in any two given directions is same, qkpz- Thus, 2A^ pairs 
characterized by qkpz are required to define the morphol- 
ogy of a surface whose growth dynamics is governed by 
KPZ equation on a three dimensional surface. Consider 
the linear couplings defined over this three dimensional 
surface. There are 3A'' linear couplings on the surface. As 
has been argued in the case of 2-1-1- dimensions, in the 
present case, on this KPZ surface only 2A'' cross couplings 
are characterized by qkpz - Both the couplings are defined 
over the set of sextet of steps defined at A^ sites. Follow- 
ing the arguments as in the case of (2-1-1) dimensional 
KPZ growth one arrives at the relation Skpz = ^Qkpz on a 
three dimensional surface, s = {2/3)q gives a — 0.28125. 

By continuing the argument for 2-|-l- and 3-1-1- dimen- 
sional growth it is possible to obtain the relation be- 
tween Skpz and qkpz in any dimension. On a d dimen- 
sional surface, a can be determined from the condition 
Skpz = ^^Ikpz which provides next higher value of a. 
Thus for KPZ equation, a can be determined in any di- 
mension. This implies that there is no upper critical 
dimension for the KPZ equation. For EW term on a sur- 
face with dimension d > \, the condition is always s = q 
giving a = 0. 

Next we consider the application to fourth order terms. 
It has been shown [16] that the relevant non linearities to 
this order are V^(V/i)^ and V • (V/i)^. The correspond- 
ing linear term is V^ft., referred as MuUin's term [17]. 
The application of perturbation up to second order to 
V^(V/i)^ term shows that its coefficient is not renormal- 
ized [18]. In the case of V • (V/i)^ it is shown to generate 
EW term [19,20]. 

In 1-1-1- dimensions a is unity for fourth order linear 
and V^(V/i)^ term [16,18]. In the later case it implies 
that the corresponding coefficient in the growth equation 
is not renoruialized. It is 0.5 for V • (V/i)"^ term [19]. 
Consider the quantity < {5xi-i — 5xiY = 2 — 2s. It 
represents average local second derivative on the surface 
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and the difference in the local particle current for the EW 
term. If this difference is not zero on the average, then 
corresponding term contributes in the growth evolution 
of the surface. In the present case it means that EW term 
is present if this difference is not zero. Hence for any 
other conservative term to be exclusively present in the 
growth equation, this current difference must be zero. It 
is zero for s = 1 giving a = 1. Thus in 1+1 dimensions for 
V2(V/i)2 and V^/i terms s = 1 [2,3]. The term V • {Vhf 
is similar to EW term in that the current for this non 
linear term is related to EW current and preserves the 
sign of the current locally. Thus condition for maximum 
number of configurations for EW term is same as that for 
the non linear term. Thus for the V- (Vh)^ term s = or 
a = 0.5 [19] in 1+1 dimensions. Thus in 1+1 dimensions 
with non conserved noise, s = and s = 1 are the only 
two possibilities defining the surface morphology. 

In 2+1 dimensions as before we must consider linear 
and other than linear (oti) couplings. So far we have 
encountered cross coupling as one of the otl couplings. 
The main difference in the analysis of the higher order 
terms using these couplings compared to second order 
terms is that, 1) more than one lattice points are to be 
considered for defining the discrete form of the term and 
2) otl couplings other than cross couplings arc needed. 
The roughness exponent is obtained by equating the total 
minimum linear couplings required to define the local 
morphology to the total minimum otl couplings along 
with the Eq. (12). 

Consider the term V^(V/i)^ or Idv term, in 2+1- di- 
mensions. In this case the corresponding discrete form 
is ^'^Xi-^-lJ + d^Xi-i,j — AS'^Xi^j + S'^Xij^i + S'^Xij^i + 
5'^yij+i+6'^yij-i-4:5^yij+S^yi+i.j+5^yi-i.j. There arc 
ten steps required to define the local morphology of the 
surface. As can be seen from Fig. 2, these steps are spread 
over total five lattice points. The morphology at every 
lattice point must be defined over five lattice points as- 
sociated with the given lattice point. To define the local 
morphology in terms of linear couplings, we note that lin- 
ear couplings at (i-l,j), (i,j+l), (i+l,.i) and (i,j-l) arc nec- 
essary. These contain total of 16 steps. However, the very 
nature of linear couplings is such that minimum two pairs 
of steps in orthogonal directions are needed at a given lat- 
tice point to allocate it the correct relative height. Thus 
the number of steps involved are in excess to the number 
required to define the local morphology. This situation is 
also encountered in the case of KPZ equation previously. 
The contribution of linear couplings is 8s. Let's consider 
otl couplings. Since the linear couplings are added, in 
the explicit form, all the steps, defining the morphology 
appear in the expression involving sum of c;ight different 
linear couplings in terms of the steps. Hence, while otl 
couplings are considered, the relative sign of all the cou- 
plings therein must be same. This ensures the presence 
of all the steps in the explicit expression for otl couplings. 
This means that otl couplings be formed between steps 



that are parallel or following. The example of a following 
type coupling is < dyi-ijSxij > where the coupling is 
formed between steps that are following (arrows) as in 
Fig. 2. This coupling is qS"^. To discriminate from the 
other cross coupling such as < Syij6xij > which is —q6^, 
we will call this coupling proportional to — g as dispers- 
ing one. We have used this later type in arriving at the 
results for EW and KPZ equations. In the present case 
following otl couplings are considered. < Syi-ijSyij > 
■ < Syi+i,jSyi^j >, < Sxij^iSxij >, < Sx^j+iSxij >, 
< 6xi-ijdyij+i >, < dyij-idxi+ij >. There is an- 
other possibility of addition of five dispersing cross cou- 
plings from Fig. 2. These are < (Jj/j+i jJxi+i.j >, < 
5xij-i5yij-i >, < Sxij+iSyij+i >, < 6xi-ijSyi-ij > 
and ,< 6xij6yij >. These two possibilities are chosen 
because they are additive involving all the steps needed 
to define the V'^iyhY term. The couplings also are dis- 
tributed to render correct weightages to the lattice points 
involved. Thus from Fig. 2 it is seen that points at 
(i+l,j) and (i,j+l) contribute three steps each while (i- 
l,j) and (i,j-l) contribute two each. Thus in both the 
sets, ratio of the couplings at (i+l,j) and (i,j+l) to those 
at (i-l,j) and (i,j-l) is 3:2. Any other additive possibilities 
do not render this ratio. It is further required that in the 
otl couplings height at (i,j) must appear. This require- 
ment is that of locality of the term V^(V/i)^. Consider 
e.g. < 6xij+i6yij+i >. This coupling does not contain 
h(i,j). Such a coupling will independently couple with 
the steps at (i,j+l) leading to second order term. We 
will illustrate this point in connection with fourth order 
linear term. This requirement eliminates set of five dis- 
persive couplings. As can be verified from the Fig. 2, all 
the couplings in other set involve h{i,j). 

In order to determine surface corresponding to the otl 
couplings preserving the locality, we need to find parallel 
coupling of the form < 6xij-i6xij >. This is easily 
obtained by noting that = {h{i, j) — h{i + + h{i + 
l,j)-hii + l,j-l) + h{i + l,j-l)-h{i,j-l) + h{i,j- 
1) — h{i,j)). Squaring this expression and writting in 
terms of steps gives the desired result, 

< Sxij-iSxij >= (1 + 2q)5'^ 

We denote this coupling as qp. < 6xi-ij6yij+i > can 
be obtained by noting that 5" = {6xi-ij + 6yij + 
6yij+i)^/6^. This leads to 

< Sxi-ijSyij+i >= (5" - 3 - 2g - 2s)/2) 

We denote this coupling as g,,. We note that \/5 = 
2i. 160964 _ rpj^g relation between the linear and otl cou- 
plings is obtained from, 

8s = Aqp + 2q, (14) 

Using Eq. (11), (12) and (14) we obtain q = —0.206 giving 
a = 2/3 with possible numerical error beyond 4th digit. 
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Thus for the Idv term present analysis shows that there 
are two possible surfaces, one with a = 0.303 and the 
other with a = 2/3. The scaling transformation gives 
a = 2/3 [16]. This shows that the system has a fixed 
point where a = 2/3 which also happens to be the value 
obtained from the scaling transformations. Under these 
conditions, since noise is not rcnormalizcd imdcr con- 
served growth conditions [3], the coefficient of the Idv 
term does not get renormalized. This shows clearly why 
perturbation analysis [18] does not show rcnormalization 
even up to second order. Thus for Idv term the exact 
value of a and the one obtained from scaling transforma- 
tions coincide. The situation in this respect is similar to 
any linear term. 

To obtain the roughness constant for V^/i refer to Fig. 
3. All the necessary steps contributing to the discrete 
form of this term are included in the diagram. By inspec- 
tion of the Fig. 3, there are linear couplings at (i-|-l,j), 
(i-l,j), (i,j+l) and (i,j-l) as before contributing total of 
8s to the linear coupling. It is interesting to verify the 
condition of locality and effect of its non observance in 
this case. From the Fig. 3, it is easily verified that one 
can obtain exactly eight cross couplings with six of them 
not involving the height h{i,j). We immediately obtain 
the condition s = q which is same as second order EW 
term. Thus it shows that the cross couplings indepen- 
dently get coupled to different lattice points affecting 
heights at those points independently. This results in 
to effectively a single lattice point after discretization. 
However, single lattice point description corresponds to 
second order term. Hence, the coupling of h{i,j) is nec- 
essary in otl couplings. Using the valid otl couplings, 

8s = Sqp + (15) 

gives 5 = or a = 1 for the unique surface corresponding 
to the term V^/i [2]. 

Finally we consider the term V • (V/i)^. From the scal- 
ing transformations, in 2+1 dimensions, a = 1/2 for 
theis term. The steps included in its discrete form are 
shown in Fig. 4. The linear couplings contributing to 
the morphology arc at (i,j) and (i+l,j+l) lattice points. 
Thus the contribution from linear couplings is 4s. There 
are two possible sets of otl couplings in this case satis- 
fying the condition of additive couplings, proper lattice 
point weightages and locality. First one will consist of 

< 5xi-ij6yij>, < 6xijSyij-i >, < Sxij5xij+i >, and 

< Syij5yi^i_j >. This will contribute 4s = 2q + 2qp. 
This results in to g = 0.25 or a = 0.58497. The other set 
is < 6xi-ij5yij >, < 6xijSyij-i >,< Sxijdyi^ij >, 

< Sxi_j^iSyij >. This will contribute 4s = Aq thus 
leading to EW term [18,20]. Our analysis shows that 
this term has additional fixed point corresponding to 
a = 0.58497. The simulations and perturbation approach 
confirms to only EW behaviour [18,19]. Significance of 
this additional fixed points needs further attention. 



We will now consider all the non linearities derivable 
from the non Inear terms considered above. Consider 
terms of the form V • (V/i)^"+^ with n=l, 2, 3, ... . The 
relevant diagrams representing these terms in the discrete 
form are exactly same as that corresponding to the term 
with n=l. Hence, all such non linearities will have same 
behaviour i.e. like EW term [20] and with additional 
surface possible having a = 0.58497. Consider terms 
V^(V/i)^" with n= 1, 2, 3, ... . In this case also, we find 
that the diagrams representing the discrete form arc all 
same as that corresponding to n=l. Thus the behaviour 
of all such non linear terms is same as that correspond- 
ing to the Idv term. For non conservative terms such as 
(V/i)^" with n=l, 2, 3, ... we find that all such terms 
will behave as the KPZ term asymptotically. These re- 
sults show that in the conserved or non conserved growth, 
the EW, KPZ and Idv are the only relevant terms up to 
fourth order. Further non linearities if relevant must con- 
tain higher number of derivatives than in Idv term STich 
as V^(V/i)'^ etc.. Experimental values of exponents other 
than the terms discussed here are likely to be the result 
of correlated noise or transients in the growth. 

In conclusion, we have shown that using the discrete 
form of the terms in a growth equation and the use of 
linear and cross coupling between steps at a site, a can 
be exactly determined for various terms in 1-1-1- dimen- 
sions, and in 2-|-l- dimensions. In particular for KPZ 
term, the existence of a in all the higher dimensions 
shows that there is no upper critical dimension for the 
KPZ equation. The method can be extended to any non 
linear stochastic equation with non conserved noise. Ex- 
act results for fourth order non linearities are obtained 
and show that the non rcnormalization for Idv term is 
a consequence of the coincidence of the exact value of 
a with that obtained from the scaling transformations. 
This study indicates that in the experiments the results 
arc likely to be affected by the presence of the difFc^rcmt 
non linear terms , the correlated noise and the transients. 
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FIG. 1. Slope distributions corresponding to s = (dotted 
curve) and s = 1 (solid lines). 
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FIG. 2. Steps contributing to the term V^(V/i)^. Each 
step is represented by an arrow as per the definition given in 
the text. 
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FIG. 3. Steps contributing to the term V*h 
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FIG. 4. Steps contributing to the term V • (V/i)^ 
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